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QCD with chiral 4-fermion interactions (xQCD)* 
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Lattice QCD with staggered quarks is augmented by the addition of a chiral 4-fermion interaction. The Dirac 
operator is now non-singular at niq = 0, decreasing the computing requirements for light quark simulations by at 
least an order of magnitude. We present preliminary results from simulations at finite and zero temperatures for 
TUq = 0, with and without gauge fields. 



1. INTRODUCTION 

xQCD is QCD with a chiral 4-fermion interac- 
tion (Nambu-Jona-Lasinio — Gross-Neveu) . 
Since the 4-fermion term is an irrelevant opera- 
tor, this theory should lie in the same universality 
class as regular QCD. 

On the Lattice we use staggered quarks. For 
simplicity we consider a theory where the 4- 
fcrmion operator has the 

C/(l) X U{1) c SU{Nf) X SU{Nf) 

flavour symmetry generated by (1, 175^5) j^]. The 
molecular dynamics Lagrangian for the lattice 
theory is 

L = -/3^[l-iRe(TrnJ7C/C/C/)]+^V'UUV' 
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where 
A=]/) 



16 ^ 



(2) 



with e = (—1)^+?/+^+*. This describes 8 flavours. 
For Nf which is not a multiple of 8 we use "noisy" 

*Talk presented by D. K. Sinclair at LATTICE'96, 
St. Louis MO, 4-8 June, 1996. 

tSupported in part by NSF grant NSF-PHY92-00148. 
tSupported by DOE contract W-31-109-ENG-38. 



fermions and multiply the fermion kinetic term 
by iV//8. 

The Dirac operator of regular lattice QCD be- 
comes singular as 0, whereas that for 
xQCD remains non-singular at rUq = 0. Conju- 
gate gradient inversion of the regular QCD Dirac 
operator requires a number of iterations which 
— > 00 as vriq — > 0. Inversion of xQCD Dirac op- 
erator requires a finite number of iterations even 
at niq = 0. In practical terms, conjugate gradi- 
ent inversion of the regular QCD Dirac operator 
requires ^10,000 iterations at the physical u and 
d quark masses. Simulating xQCD at niq = 0, 
the worst we have found is ^ 400 iterations — a 
saving of a factor of ~ 12 in CPU time! xQCD al- 
lows us to work at the chiral limit, which is often 
useful. 

All the simulations which I will now describe 
are being performed at mq — 0. Simulations are 
being performed on the CRAY C-90 and CRAY 
J-90 at NERSC. 



2. xQCD AT FINITE TEMPERATURE 

We are studying the finite temperature be- 
haviour of xQCD with Nf — 2 and rUq — 0, 
using the hybrid molecular dynamics algorithm 
with "noisy" quarks, as a function of /3 and 7. 

We are studying the deconfinement and chiral 
transitions — not in general coincident at finite 
a. The 2 transitions come together as a ^ and 
in general as 7 is increased. By measuring the 
critical behaviour of such observables as (ipip) , we 
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hope to determine the equation of state, where 
these transitions coincide. For this, being able to 
work at the chiral Hmit is a boon. 

Our present simulations are being performed at 
Nt = 4, using 8'^ X 4 lattices to find the position 
of the phase transitions, and 12^ x 24 x 4 lattices 
to study the transition in detail (the extension 
of this lattice in the z direction is to allow the 
measurement of hadron screening masses). 

At 7 = 2.5 the deconfinement transition occurs 
between /3 — 5.5 and /3 — 5.6 while the chiral 
transition is at /3 ~ 10. At 7 = 5 the transi- 
tions still appear separated, but clearly the chiral 
condensate feels the deconfinement transition. By 
7 = 10 the 2 transitions are close enough together 
that it is no longer clear whether they are distinct 
- see figure |[ The way we measure (iptp) , it does 
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Figure 1. Wilson line and (ipip) as functions of p 
for 7 = 10 

not vanish in the symmetric phase, but should 
instead behave as 1/ \/V. We have tested this at 
'Y — 10, /3 — 5.4 and (3 — 5.35, and find it to 
be true within errors which indicates that both 
these /3's lie in the chirally symmetric ~ decon- 
fined - phase. (The phase transition occurs at 
(3 « 5.33). 



3. ZERO TEMPERATURE xQCD 

We are simulating xQCD with Nj = 2 
ruq = at 7 = 10 and (3 — 5.4 on an 8^ 



and 
X 24 



lattice to study zero temperature physics and in 
particular hadron spectroscopy. 

We have measured the propagators of the a 
and TT auxiliary fields on 2500 configurations sep- 
arated by 2 time units. The first 500 were dis- 
carded leaving 2000 for analysis. These were 
binned into 200 bins of 20 time units in length 
for error analysis and fitting. The pion propaga- 
tor was fitted to the massless scalar lattice prop- 
agator (see figure H). The fit, over time interval 
1-12, had confidence level 78% showing that it 
is indeed a Goldstone pion. The a propagator is 
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Figure 2. The propagator for the Goldstone pion. 
The fit is to a massless scalar lattice propagator. 

noisy, as expected. The best fit we were able to 
obtain gave a mass of 1.12(23). Current statistics 
preclude determining whether there is any signif- 
icant contribution from the 2-pion cut. 

So far we have stored 250 configurations spaced 
by 20 time units for calculating the hadron spec- 
trum. We plan first to calculate the local meson 
and baryon spectra. For the mesons, this will in- 
clude the normal, chirality 1 mesons (quarks hav- 
ing chirality ±i) formed from linear combinations 
of 



Tr[G'(7r)Gt( 



(3) 
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and the chirality mesons formed from linear 
combinations of 

Tr[G(7r)Gt(^)] (4) 

where G is the quark propagator. For the nor- 
mal meson propagators, the a and tt have dis- 
connected contributions which we will also need 
to calculate, while the p and ai do not. For the 
chirality mesons there are no disconnected con- 
tributions. 

Working at ruq = we will need a careful finite 
size analysis of the hadron spectrum, since hadron 
sizes are determined by m^r. 

4. xQCD AT ZERO GAUGE COUPLING 

The /3 ^ oo limit of xQCD with Nf = 
2 is the 4-d Nambu-Jona-Lasinio — Gross-Neveu 
model with Nf = 6. We have studied this the- 
ory on an 8"* lattice to find its chiral transition, 
which enables us to identify its strong and weak 
coupling domains. We find the transition 7 to be 
~ 1.7. We have also studied its finite temper- 
ature behaviour on an 8^ x 4 lattice, where the 
transition occurs at 7 ^ 1.5. 

5. COMMENTS ON xQCD AT FINITE 
BARYON NUMBER DENSITY 

Quenched lattice QCD has problems at finite 
baryon number density. If /i is the quark num- 
ber chemical potential, one would expect QCD to 
have a phase transition at /i w toat/B. Instead, at 
rUq = 0, quenched QCD appears to have a tran- 
sition at ^ = 1^. 

The nearest thing to a quenched version of 
massless xQCD we know is quenched QCD with 
fermions interacting with a mean field approxi- 
mation to the chiral fields in which a is replaced 
with a constant and the tt with zero. This is just 
massive quenched QCD (mass =< a >) at fi- 
nite chemical potential, which has a transition at 
jj, w m7r(< cr >)/2 which is no longer zero. Since 
we expect that m.,r{< cr >)/2 < mtq{< a >)/3 
we have not completely solved the problem. How- 
ever, for moderate < cr >, the behaviour is more 
physical than for conventional quenched QCD 
— the unphysical region from mT(< a >)/2 to 
r7ijv(< cr >)/3 is small. 



6. SUMMARY AND CONCLUSIONS 

xQCD enables simulations at physical u and d 
quark masses with at least an order of magnitude 
saving in CPU time. It also enables simulations 
with zero quark masses which is important for 
determining the equation of state. A renormal- 
ization group analysis will be needed to continue 
to the continuum limit. 

After finishing this first round of simulations 
and hadron mass/screening mass measurements 
at zero and finite temperature, we will move to 
larger lattices, and also perform some simula- 
tions at non-zero quark mass. The next stage 
will involve using a more physical chiral group 
iSU{2) X S'[/(2)??). 

Auxiliary fields and constituent quark masses 
give promise of an effective Lagrangian interpre- 
tation of xQCD. These effective quark masses 
could (partially) resolve the problems at finite 
baryon number density. 

Similar theories have been considered by 
Brower, Orginos, Shen and Tan Theories with 
4-fermion and gauge interactions have a long his- 
tory; see for example 0. 
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